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We theoretically study equilibrium and dynamic properties of nanosized magnetic skyrmions in 
thin magnetic films with broken inversion symmetry, where electric field couples to magnetization 
via spin-orbit coupling. Based on a symmetry-based phenomenology and micromagnetic simulations 
we show that this electric-field coupling, via renormalizing the micromagnetic energy, modifies the 
equilibrium properties of the skyrmion. This change, in turn, results in a significant alteration of 
the current-induced skyrmion motion. Particularly, speed and direction of the skyrmion can be 
manipulated by designing a desired energy landscape electrically, which we describe within Thiele’s 
analytical model and demonstrate in micromagnetic simulations including electric-field-controlled 
magnetic anisotropy. We additionally use this electric-field control to construct gates for controlling 
skyrmion motion exhibiting a transistor-like and multiplexer-like function. The proposed electric- 
field effect can thus provide a low energy electrical knob to extend the reach of information processing 
with skyrmions. 
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I. INTRODUCTION 

One of the central challenges faced by the field of spin- 
tronics, as an alternate to conventional charge-based in¬ 
formation processing, is to encode information in a mag¬ 
netic configuration (i.e. a collection of spins) which is 
both stable against thermal fluctuations and electrically 
manipulated with minimal energy budget. Recently, 
magnetic skyrmions have emerged as one such candi¬ 
date configuration, opening a new sub-field of informa¬ 
tion processing with skyrmions, dubbed as skyrmionics 
Skyrmions are topologically nontrivial magnetic textures 
characterized by a non-zero integer-valued skyrmion 
number A/*, which counts the number of times the tex¬ 
ture wraps a sphere in the spin-spac^. Much recent in¬ 
terest has been fuelled by observation of the skyrmion 
crystal phas^, its emergent electrodynamic propertied 
and manipulation by currents at low temperatures . 5 

Meanwhile, motivated by room temperature ap¬ 
plications, thin films with perpendicular anisotropy 
have been proposed as hosts for controlled nucle- 
ation and manipulation of individual skyrmions 5 . Sin¬ 
gle skyrmions in such films can be stabilized by ei¬ 
ther Dzyaloshinskii-Moriya interaction(DMI^ or dipole- 
dipole interaction(DDI)Competing against Zeeman and 
exchange energies. Additionally, these films can be 
typically sandwiched between a dielectric and a high 
spin-orbit coupling material resulting in broken inver¬ 
sion symmetry along the growth direction. Such film 
structures, apart from allowing for a finite DMI, have 
recently been shown to exhibit current-induced spin- 
orbit torques(SOT^P, which, in turn, result in efficient 
current-induced motion of chiral magnetic texturesP 
These attractive properties have made such films an 
ideal candidate for skyrmion-based race-track-like mem¬ 
ory applications . 5 9 In fact more recently, experimental 
evidence for skyrmionic configurations stabilized at room 
temperature in thin films 10 and multilayers 11 in proxim¬ 
ity with heavy metals have been presented. In Ref. HI 


current-induced motion has also been observed. 

In this article, we propose and show—via a combina¬ 
tion of symmetry based phenomenology and micromag¬ 
netic simulations—that an electric-field applied across 
such film structures can dramatically influence current- 
induced skyrmion motion. In particular, we find that 
for a fixed current, the electric-field can control both the 
speed and the direction of skyrmion motion. More im¬ 
portantly, we show that by applying electric-fields in a 
certain pattern, defined via gated structures, skyrmions 
can be guided along desired trajectories. We also identify 
the mechanism for this electric-field control of skyrmion 
motion, which can be explained in terms of the change in 
the static properties of skyrmions, brought about by the 
electric-field induced modification of the micromagnetic 
free energy. Finally, utilizing this additional electrical 
knob we show that the realm of skyrmionics can be ex¬ 
tended beyond the race-track-like memory applications, 
to possibly realize logic and computing functionalities. 


II. ELECTRIC-FIELD COUPLING: 

PHENOMENOLOGY 

We begin by constructing a general phenomenological 
model describing electric-field coupling to the magnetiza¬ 
tion dynamics based on the symmetries of the structure. 
We are interested here in two dimensional ferromagnetic 
films considered to have the following symmetry proper¬ 
ties (see Fig. fTI) : broken inversion symmetry along the 
growth axis (oriented along z ), by interfacing it with 
different materials on both sides and application of an 
electric field, which is assumed to be applied along the 2 
axis via gates. While for simplicity, the structure is con¬ 
sidered isotropic about the z axiff 2 . In addition, we are 
interested in the temperature regime well below the Curie 
temperature of the ferromagnet, where, suppressing the 
magnitude fluctuations of magnetization, the magneti¬ 
zation dynamics can be expressed within the Landau- 
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FIG. 1. Schematic of the system: A ferromagnet(FM) layer 
with a skyrmion, sandwiched between a heavy metal(HM) 
and a dielectric, breaking the symmetry along the z axis. A 
current through HM with density J drives skyrmion via SOT, 
while simultaneously local electric fields (oriented along z) can 
be applied at a gate region by an external voltage V . 


Lifshitz-Gilbert framework as 13 : 

<9/111 = 7in x H* + am x < 9 t m. (1) 

Here m is a unit vector along the magnetization M 
(with magnitude M s ), i.e. M = M s m, while 7 and a 
label gyromagnetic ratio and Gilbert damping param¬ 
eters, respectively. Moreover, we have defined H* = 
—+ H( J, E), with representing the effective 

field contribution derivable from a micromagnetic free 
energy J 7 , and H(J,E ) denoting the contribution due 
to application of current (having a current density, J ) 
and electric-field(E). In its minimalistic form, the micro- 
magnetic free energy is comprised of the classical dipole- 
dipole interaction Td, the exchange interaction parame¬ 
terized by an exchange constant A, and the Zeeman en¬ 
ergy due to an applied magnetic field H a . In addition, 
the breaking of inversion symmetry about the z axis in 
the presence of spin-orbit interaction gives rise to interfa¬ 
cial perpendicular anisotropy and DMI^, whose strengths 
are parameterized by K and 19, respectively, resulting in: 


The terms proportional to t\f> and (d represent the dissi¬ 
pative components, i.e. which cannot be derived from a 
free energy, while terms proportional to t]f and repre¬ 
sent current and electric-field-induced modification of T. 
In particular, the terms proportional to current density 
represent the field-like and anti-damping like components 
of spin-orbit fields, microscopically derived from spin- 
Hall and Rashba/Edelstein effects^. Additional current- 
dependent terms representing corrections due to angle- 
dependence of spin-orbit field and symmetry-allowed dis¬ 
sipative torques inexpressible in the simple anti-damping 
like form- 5 are not included for simplicity. On the 
other hand, the parameter describing coupling of electric- 
field to magnetization, denoted by fr, can be identi¬ 
fied as electric-field dependent perpendicular magnetic 
anisotropy. Microscopically, electric-charge-^, strain- 17 
and ionic movement-inducecf 18 anisotropy modification 
are few examples which are gaining increased attention 
in the spintronics community for exploring, both, fun¬ 
damental physics and potential low-power technological 
applications 19 . Particularly for domain-wall based de¬ 
vices, the electric-field dependent anisotropy has been 
shown to trap domain walld^ and proposed to induce 
domain-wall motion 2 0. Finally, the dissipative coupling 
of electric-field should become relevant for very thin 
dielectric^l which, to the best of our knowledge, has 
not yet been measured experimentally. 

At this point we mention for completeness that terms, 
allowed by symmetry, where electric-field and current 
couple to gradients of magnetization, which could them¬ 
selves be interesting to explore in proposed experimental 
structures, have been omitted from Eq. [3] For exam¬ 
ple, the electric-field could generate a term of the form 
~ Em z V • M, describing microscopic electric-field de¬ 
pendence of DMI. Such electric-field-induced DMI has 
been proposed to induce a transverse ME effect in spiral 
spin magnets 2 ^! Motivated by including minimal ingre¬ 
dients to stabilize skyrmions and study non-trivial effects 
of electric-field on their dynamics based on the available 
experimental data, we set D = rjp = (d = 0 in the 
following sections. However, the qualitative features of 
electric-field control and its mechanism are expected to 
be extended for the case when some of these parameters 
are included as well. 


E = E d + A[{d x m) 2 + (^m) 2 ] — M H a (2) 
—Km 2 z + Dm z \7 • m. 


III. ELECTRIC-FIELD EFFECT: STATICS 


Here, V = d x x + d y y. Similarly, the current and elec¬ 
tric field-induced spin-orbit fields, obtained upon impos¬ 
ing the above mentioned structural symmetries on Eq. [T| 
(noting J and m z are even, while E , m x and m y are odd 
under inversion about z axis) and, in the spirit of linear 
response, retaining terms linear in E and J, are given by: 

H( J, E) = tjdJ m xy + (dEi71 z m x z (3) 
+VFJy + CfEuIzZ. 


The dynamic control of skyrmions, as will be explained 
later, can be understood in terms of the modification of 
their static properties. To this end, we discuss first the 
electric-field effect on the static properties of a skyrmion. 
In general, two separate regimes can be defined based on 
the stability of the skyrmion. In the first regime, the 
skyrmion remains stable with a smooth change in its di¬ 
ameter D s , and the domain-wall width A w , surrounding 
the skyrmion. This change occurs in order to minimize 
the micromagnetic energy corresponding to the electric- 
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FIG. 2. Static properties of skyrmion as a function of perpen¬ 
dicular anisotropy K : broken vertical line marks approximate 
boundary between two regions. For K > 4.25 X 10 6 erg/cm 3 , 
skyrmion is stable, with its micromagnetic configuration for 
K = 4.5 x 10 6 erg/cm 3 shown at top right corner. Here m z de¬ 
notes the out of plane component of the magnetization vector, 
while the in-plane components are shown with an arrow. The 
diameter and domain-wall width surrounding the skyrmion 
in this region is represented by triangle and square markers, 
respectively. For K < 4.25 x 10 6 erg/cm 3 the skyrmion is 
unstable, instead a multi-domain/ in-plane state is favored. 


field induced modification of perpendicular anisotropy: 
the energy cost, due to deviation of magnetization from 
the easy axis in the wall region, increases with increasing 
K. In an attempt to reduce this energy cost, the domain 
wall area (~ 2ttD s A w ) is reduced resulting in decrease 
of both diameter of skyrmion and the wall-width. In 
the second regime, on the other hand, the skyrmion it¬ 
self becomes unstable. This instability occurs close to 
the out-of-plane to in-plane transition of magnetization 
when the perpendicular anisotropy is lowered to a value, 
such that, the perpendicular anisotropy energy is over¬ 
whelmed by the dipole-dipole interaction, due to the out 
of plane component of the magnetization^. 

To demonstrate the qualitative picture presented 
above in a specific numerical model, we perform micro- 
magnetic simulations for a DDI-stabilized skyrmion(i.e. 
with DMI strength set to zero) with saturation mag¬ 
netization, M s = 920 emu/cm 3 and exchange stiffness 
Ami /ierg/cm. In order to take electric-field-induced 
modification of perpendicular anisotropy into account, 
the strength of perpendicular anisotropy was varied in 
the range K = [2-5.5] x 10 6 erg/cm 3 . This range cor¬ 
responds to films with K = 4.25 x 10 6 erg/cm 3 in 
the absence of electric-field, and a maximum value of 
electric-field being ^0.1 MV/cm, for a coupling con¬ 
stant ( = 20 erg/(V-cm 2 ^. These material parame¬ 
ters are typical for strain-based ferromagnet/dielectric 
heterostructures such as CoPd(alloys)/PMN-PT, which 


are chosen as they exhibit one of the largest values of 
magneto-electric coupling (i.e. () at room temperaturd^. 
In addition, a DDI skyrmion state has been found to 
be stable in finite-sized nanodisks of cobalt- and iron- 
based ferromagnets with similar material parameters 27 . 
The simulated magnetic films are assumed to be two- 
dimensional with periodic boundary conditions in the 
film plane and a thickness of t = 32 nm. Additionally, 
an external magnetic field of strength H a /Hd = 0.3 with 
Hd = 47rM s , is applied along the 2 axis, which is needed 
to stabilize a DDI skyrmion for films of infinite extent 2 ^. 
All micromagnetic simulations are performed using LLG 
Micromagnetic SimulatoP^ with temperature set to zero. 

To study the static properties of a skyrmion, the ini¬ 
tial magnetization was set close to a film with a sin¬ 
gle skyrmion state, i.e. with a central region of suffi¬ 
ciently large diameter pointing along —z and the rest 
of the film oriented along +z. This initial configura¬ 
tion was then allowed to relax to equilibrium for each 
value of perpendicular anisotropy. The result of such 
a procedure is summarized in Fig. [2] For anisotropy 
strength K < K c = 4.25 x 10 6 erg/cm 3 , the single 
skyrmion state is unstable, favoring a multi-domain (for 
relatively strong iF, i.e 3 < K < 4.25 x 10 6 erg/cm 3 ) 
or in-plane magnetization (for still weaker anisotropies, 
i.e. K < 3 x 10 6 erg/cm 3 ). Such a behavior is typical 
for films, where the out-of-plane magnetization-induced 
dipole energy (^ M 2 ) starts dominating the perpendic¬ 
ular anisotropy energy K). This picture can be fur¬ 
ther confirmed by calculating the so-called quality factor 
Q = K/2irMg which has a value of Q c ~ 0.8, i.e. close 
to one, at the critical anisotropy, K c . On the other hand, 
for K > K c in the simulations, the resulting equilibrium 
configuration is that of a single skyrmion (shown in the 
inset of Fig. [2| 29 . As expected, in order to minimize 
the dipole-dipole interaction energy in the wall, the DDI 
skyrmion is “Bloch-like” with skyrmion number J\f = 1. 
Moreover, increasing anisotropy results in decreasing di¬ 
ameter and wall-width surrounding the skyrmion, sup¬ 
porting the qualitative argument presented above. 


IV. ELECTRIC-FIELD EFFECT: DYNAMICS 

In this section we show how the modification of static 
properties of the skyrmion results in the electric-field con¬ 
trol of current-induced skyrmion motion. Correspond¬ 
ing to the two regimes, defined for modification of static 
properties, there exist two qualitatively different regimes 
for the dynamic control, (i) When a uniform electric field 
with a magnitude below a critical value is applied, such 
that the skyrmion remains stable, the change in D s and 
A w results in change in the magnitude and the direction 
of skyrmion motion, (ii) On the other hand, defining 
unstable regions for skyrmions by application of a non- 
uniform electric-field above the critical value, results in 
constraining skyrmion motion along desired paths (which 
we refer to as the “guiding regime” below). To gain more 
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insight into this connection between electric-field-induced 
static and dynamic control, we begin with an analyti¬ 
cal model of current-induced skyrmion motion based on 
Thiele’s collective coordinate approach!^. 



anisotropy (10 6 erg/cc) 


FIG. 3. Velocity modulation for uniform electric-field: the 
schematic in the bottom shows the steady-state of skyrmion 
moving with velocity v and the corresponding forces acting on 
it. F s represents current-induced spin-orbit force, while 0v 
and Tv represent the Lorentz and damping force, respectively. 
Open square and triangle markers represent speed and the 
angle 'ip (between current and v as shown in the schematic) 
as obtained from analytical model. The corresponding speed 
and angle from micromagnetic simulations are represented by 
the closed square and triangle markers. 

Analytical model — In the collective coordinate ap¬ 
proach, the relevant degrees of freedom, describing the 
equation of motion of a magnetic-texture, are given by 
the so- called soft modes 30 . The most relevant soft modes 
for a skyrmion in an infinite film are its rigid translations 
in the film plane. Consequently, the collective coordi¬ 
nates included for describing the motion of a skyrmion 
are (X, V), which label the position of the skyrmion’s 
center along the x and y axis, respectively. Within this 
approximation, LLG Eq. fl]can be reduced to the follow¬ 
ing equation of motiomE^F 

G x v - Tv + F s - F u = 0, (4) 

where v = (X,T) is the velocity vector, while G = 
(0,0,0) with Q = f m • ( d x m x d y m). Inclusion of 
other modes, such as wall fluctuations, results in de¬ 
viations from the rigid motion approximation and en¬ 
dows Eq. [4] with a “mass term’®, which is not the fo¬ 
cus of the current study. The first term describes the 
“Lorentz force”^ 2 due to the skyrmion’s non-trivial topol¬ 
ogy induced fictitious magnetic field 33 . Furthermore, the 
integrand in Q is the local solid angle, making Q inde¬ 
pendent of the exact magnetization texture and simply 
proportional to the skyrmion number, i.e. Q = 4 irAf. 


The second and third terms originate from the Gilbert 
damping and current-induced spin-orbit field term in 
LLG Eq. [I] respectively, with T = a f d x m • d x m and 
Ff = 777 Jj dim • (m x y), where i labels the cartesian 
coordinates. Finally, the fourth term represents the force 
due to variation of the skyrmion’s “potential energy” [/, 
i.e. F u = —VC/. The skyrmion’s potential energy is, in 
turn, defined as U(X,Y) m f T + CpEra 2 /2, where the 
integral is performed for the magnetic configuration of a 
skyrmion with its center located at (X, Y). In contrast 
to 0; T, F s and F u depend on the exact functional form 
of the magnetization. 

In general, the functional form of F u will depend 
on the interaction of the skyrmion with its “environ¬ 
ment” and will be discussed for the two regimes below, 
while analytical expressions are obtained for the dissi¬ 
pative and spin-orbit forces using Thiele’s ansatJ^ 4 ! for 
a skyrmion’s magnetic configuration as T = a7r 2 D s /A w 
and F s = ( 0 , 7 ? 7 J 7 r 2 D s /2). Within this Thiele’s ansatz, 
the integrand for T and F s is assumed to be non-zero 
near the wall region, where the magnetization profile, pa¬ 
rameterized by the polar(#) and azimuthal(</>) angles as 
m = (sin#cos</>, sin#sin 0 , cos#), is given by cos#(r) = 
tanh( 7 rr/A w ). Additionally, in order to compare to mi¬ 
cromagnetics, the Bloch-like DDI skyrmion configuration 
with M = 1 is used, i.e. = p + 7r/2 with (r, p) rep¬ 

resenting the two-dimensional film plane in cylindrical 
coordinates. 

We begin by looking at the regime when a uniform 
electric-field, having a magnitude below the critical value 
required to make the skyrmion unstable, is applied to 
the infinite film. In this case, the skyrmion’s potential 
energy is independent of its location and hence F u = 0. 
Consequently, solving Eq. [4] for velocity, we obtain: 


v = 


_ yqJP 3 _ 

8^1 + (anD s /AA w ) 2 


; ip = tan 1 (airD s /4A w ), 


( 5 ) 

where v is the skyrmion’s speed and is the angle be¬ 
tween the skyrmion’s motion and the current (see inset of 
Fig. 3] for a schematic). This constitutes one of the main 
results of the analytical model, relating the skyrmion’s 
static and dynamic properties. The electric-field control 
of skyrmion’s motion in this regime is thus expected to 
enter via modification of D s and A w , as discussed in 
section III, and is shown in Fig. [ 3 ] for the diameter and 
domain wall width extracted from static micromagnetic 
simulations. These analytical results are compared later 
against the velocities and angles extracted directly from 
dynamic micromagnetic simulations. 

The dependence of U on electric-field suggests another 
qualitatively different route to control the skyrmion’s dy¬ 
namics via application of non-uniform electric-field and 
consequently “engineering” desired forces F u . Inter¬ 
estingly in multiferroics, electric-field gradient-induced 
forces have been theoretically shown to induce Hall-like 
skyrmion motion 35 . Here, we propose the scenario where 
unstable regions [such as the one shown in the inset of 
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FIG. 4. Skyrmion guiding regime: (a) the bottom schematic shows the steady-state of a skyrmion moving along the boundary 
between stable and unstable regions (marked by a horizontal broken line) and corresponding forces. The angle between current 
and the normal to this boundary is represented by /3. In addition to current-induced force, the damping force and the Lorentz 
force (as also shown in Fig. [3j , there exist a dipolar force between the skyrmion and in-plane magnetization in the unstable 
region, denoted by F u . The skyrmion’s speed along the boundary as a function of [3 as obtained from the micromgnetic 
simulations is shown by the square marker. While the broken curve shows the result from the analytical model, i.e. Eq. [6] for 
same current density. The origin of the dipolar force is illustrated in (b), which compares the interaction of a skyrmion with 
the environment in the presence (top schematic) and the absence (bottom schematic) of an unstable region. The red circular 
region represents the core of the skyrmion, with the magnetization inside and outside the skyrmion pointing along - and + z, 
respectively, while the lines emanating from the cylinder represent the dipolar field due to the core region. The magnetization 
in the unstable region for the top schematic points in the plane as shown by the arrows. 


Fig. [4ja)] can be defined for skyrmions, where magnetiza¬ 
tion deviates from the out-of-plane direction, due to the 
application of an electric field above the critical value. 
In this case, F u becomes non-zero as the skyrmion ap¬ 
proaches the boundary between the stable and unstable 
regions due to dipolar interaction, as explained next. The 
dipolar fields originating from the “core” (i.e. the region 
inside the skyrmion’s diameter) point in a direction op¬ 
posite to the core’s magnetization in the region outside 
the skyrmion’s core. Deviation of magnetization from 
this dipolar field direction in the unstable region will 
thus result in a higher dipolar energy when compared 
with skyrmions interacting with the film in the absence 
of unstable regions [as illustrated in Fig. Qb)]. Conse¬ 
quently, as the skyrmion moves closer to the boundary 
between the stable and unstable regions, the dipolar en¬ 
ergy increases and results in F u preventing the skyrmion 
to enter the unstable region. Due to this additional F u 
a new steady state is expected. This situation is very 
similar to the case of a skyrmion driven in nanowires, 
where in steady state the skyrmion is driven along the 
nanowire, with the dipolar forces from the boundary of 
the nanowire balancing the Lorentz forcd^. Thus generi- 
cally it is expected that, for a range of drive currents, the 
skyrmion can be guided along the boundary between the 
stable and unstable regions, as depicted schematically in 
f^. @ The advantage here, as compared to the nanowire 
case, is that this boundary is defined by application of 


the electric-field and hence can be constructed on the 
fly, providing a controllable method to “rewire skyrmion 
current”. A general analytical expression for the velocity 
depends on the details of F u , however a simple expression 
can be derived for the case when the magnetization con¬ 
figuration in the unstable region has translational sym¬ 
metry along the boundary [such as the one shown in the 
schematic of Fig. |4ja)] . This translational symmetry de¬ 
mands that F u is oriented perpendicular to the boundary, 
whose magnitude is denoted by F u . Solving for a steady 
state with the skyrmion moving along the boundary, tak¬ 
ing the angle between the current and the normal to the 
boundary to be /?, and balancing the forces along and 
transverse to the boundary [as depicted in the inset of 
Fig. |4ja)] , we obtain for the guiding regime: 

V = cos/3; F u = Anv. (6) 

2a 

We can see from the above expression that by orient¬ 
ing the electric field gate boundary with respect to the 
current, i.e. changing /?, the skyrmion’s speed can be 
controlled dramatically. The position of the skyrmion 
with respect to the boundary, on the other hand, is 
given by the balance of Lorentz and dipolar forces, with 
the skyrmion pushed closer to the boundary for larger 
currents. We will now compare these analytical in¬ 
sights to the micromagnetic simulations and, in partic¬ 
ular, prove the existence of the electric-field-controlled 
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skyrmion guiding regime. 

Micromagnetics — The micromagnetic simulations 
were performed for a DDI skyrmion with the same ma¬ 
terial parameters as used in the case for static simula¬ 
tions. In addition a current-induced spin-orbit torque, 
with gd = hJ6 s /2eM s t and 0 S = 0.2^, was turned on 
as the driving force for all the simulations presented in 
this section.® Here h , e and J are the reduced Planck’s 
constant, charge of an electron and current density re¬ 
spectively. Experimentally, this can be achieved, for ex¬ 
ample, by interfacing the ferromagnet in the ferromag- 
net/dielectric heterostructure with a heavy metal like Ta, 
Pt, W etc S3 The dynamics were calculated by integrat¬ 
ing Eq. [I] with a damping parameter a = 0.05. In Fig. [3] 
we first compare the results for a uniform electric-field 
with a magnitude below the critical field. These results 
are obtained for a single skyrmion moving in a film, with 
periodic boundary conditions along the film plane, and 
a fixed current density of J = 10 8 A/cm 2 , while simul¬ 
taneously varying the anisotropy throughout the film in 
the stable range (as obtained in the static simulations, 
i.e. for K = [4.25 — 5.5] x 10 6 erg/cm 3 in Fig{2|. As can 
be seen, both the speed and the angle of the skyrmion 
with respect to the current direction are modulated by 
application of the electric-field. Moreover for the range 
of electric-field magnitudes studied here, the trend pre¬ 
dicted by the analytical model fits the obtained v and ip 
well, both qualitatively and quantitatively, substantiat¬ 
ing the analytical picture presented above. For the mate¬ 
rial parameters used, the electric-field modulation shown 
in this regime is, however, relatively weak. For a much 
stronger control of skyrmion dynamics, we next turn to¬ 
wards the other regime, i.e application of a non-uniform 
electric-field. 

To demonstrate the skyrmion guiding regime by defin¬ 
ing unstable regions for skyrmions, micromagnetic simu¬ 
lations were performed with the value of K = 4.25(2) x 
10 6 erg/cm 3 in the stable (unstable) regions marked in 
Fig-H The dynamics was then simulated for a skyrmion 
driven by current, with the initial state of a single 
skyrmion in the stable region. In Fig. [5] we show various 
snapshots of the results of these simulations. First note 
that, as was discussed in the section III, the magnetiza¬ 
tion in the unstable region goes to an in-plane state. In 
the initial stage of the simulation, i.e. when the skyrmion 
is relatively far from the the unstable region [Fig [5jb)], 
the skyrmion behaves very similar to F u = 0 case mov¬ 
ing at an angle with respect to the current. However, as 
the skyrmion approaches the boundary between the sta¬ 
ble and unstable region [Fig.|5jc)] a non-zero F u deflects 
the skyrmion, preventing it from entering the unstable 
region. After t ~ 40ns a new steady state is reached with 
the skyrmion moving along the boundary between the 
stable and unstable region [Fig. [5|d)] and hence proving 
the existence of the guiding regime. The robustness of 
this steady state was checked by varying the magnitude 
of drive current density. For this example, the skyrmion 
is still driven along the boundary for current density be¬ 


low J ~ 5 x 10 8 A/cm 2 . As the drive force increases, 
with increasing current density, the skyrmion is pushed 
closer to the boundary according to Eq. [6j However, 
above J ~ 5 x 10 8 A/cm 2 the drive force overcomes F u , 
resulting in the skyrmion entering the unstable region 
where it is annihilated. Motivated by applications, we 
next discuss how this regime of skyrmion guiding could 
provide new opportunities for constructing possible spin- 
tronic devices. 



FIG. 5. Micromagnetic simulations showing the guiding by 
electric field: (a) Initial state showing a stable skyrmion with 
K — 4.25 x 10 6 erg/cm 3 in the region above the horizon¬ 
tal dashed line and in-plane magnetization below it, where 
K — 2 x 10 6 erg/cm 3 . A current, with density J, and 
making an angle [3 with the normal to the boundary be¬ 
tween stable and unstable regions, provides the driving force 
for the skyrmion. The inset shows the color code for out 
of plane magnetization, while the in-plane components are 
shown by arrows. The snapshot of simulation after (b) 5 ns, 
(c) 10 ns and (d) 20 ns, are shown with the arrows alongside 
the skyrmion representing the instantaneous direction of the 
skyrmion’s motion. 


V. SKYRMIONICS BEYOND RACE-TRACKS 

Nanosized skyrmions can be used as non-volatile 
information carriers in a spintronic device owing to 
their stability against thermal fluctuations. The in¬ 
formation stored in the skyrmion can, in turn, be 
read both optically or electrically (by utilizing vari¬ 
ous magneto-resistance effects such as topological Hall 
effectP^, tunneling-® or giant-® magnetoresistance). 
One such proposed example of a spintronic device is a 
race-track memory, where information stored in a mag¬ 
netic texture is moved by current to a desired location 
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where it can be read or written 41 . With the discov¬ 
ery of current-induced skyrmion motion, skyrmions in 
a nanowire are considered as a promising alternate to 
“conventional” domains in a race-track memory as they 
show very weak pinning by external impurities requiring 
orders of magnitude lower critical depinning currents 5 . 
The new possibility of controlling “skyrmion current” 
by the electrical-field presented here is thus expected to 
open up novel avenues of information processing with 
skyrmions. We demonstrate this by specifically showing 
via micromagnetic simulations that electric-field gates 
can be used to turn “skyrmion currents” off, i.e. provid¬ 
ing a transistor-like operation, and deflect information 
carried by skyrmions in specific parts of the spintronic 
circuit, i.e. a multiplexer-like operation. 



FIG. 6. Micromagnetic simulations showing the transistor¬ 
like function: top schematic shows the initial state with a 
skyrmion located on the right end of the nanowire, along with 
a current pulse of 10 ns passed through the heavy metal un¬ 
derlayer. An external electric-field can be applied locally via 
a voltage over a gate region indicated in the figure. The re¬ 
sulting location of the skyrmion at the end of the pulse when 
the gate is in OFF (bottom left) and ON (bottom right) state. 

To demonstrate the transistor-like operation, a 
nanowire geometry is adopted in the simulations with 
width of the nanowires set at 200 nm. The anisotropy 
strength was K = 4.25 x 10 6 erg/cm 3 throughout the sim¬ 
ulations everywhere in the nanowire except in a region of 
width 50 nm defined as the gate. To simulate the effect of 
electric-field applied on the gated region, the anisotropy 
value was dropped to K = 2 x 10 6 erg/cm 3 when the gate 
is in the “ON” state and kept at K = 4.25 x 10 6 erg/cm 3 
in the “OFF” state. The rest of the parameters are same 
as for the dynamical simulations presented above. The 
result of the application of a drive current pulse along the 
x axis of duration t p = 10 ns is shown in Fig. [6j When 
the gate is in the OFF state, the skyrmion is driven past 
the gate region along the nanowire (Fig. [6] bottom left). 


On the other hand, when the same current pulse is ap¬ 
plied in the presence of a gate voltage, the skyrmion is 
blocked by the gate (Fig. [6] bottom right), due to the 
repulsive dipolar forces from the unstable gated region, 
thus exhibiting a transistor-like function. This situation 
is the special case of the guiding regime, i.e. Eq. [6| with 
P = 90. 



FIG. 7. Micromagnetic simulations showing the multiplexer¬ 
like function: top schematic shows the initial state for two 
magnetic nanowires coming together at the “wing-shaped” 
junction with a skyrmion located at the right end of the upper 
arm. A local external field can be applied at the shown gate 
region. The result of the simulation at the end of a current 
pulse of 10ns passed through the heavy metal underlayer when 
the gate is OFF (bottom left schematic) and ON (bottom 
right schematic). 

Next we simulate a “wing-shaped” nanowire junction, 
with two nanowires, each of width 200 nm, meeting at 
an angle representing two different arms of a spintronic 
circuit. A gate in this case is defined on one of the arms 
as shown in Fig. [7| which also shows snapshots of the 
simulations. When the gate is OFF, a current pulse along 
the x axis in this case drives the skyrmion to the lower 
arm (Fig. [t] bottom left). This motion is similar to a 
skyrmion moving at an angle with the current due to the 
Lorentz force, as seen for a skyrmion moving in an infinite 
film. However, when the gate voltage is turned on, the 
skyrmion driven by the same current pulse is deflected off 
the gate region into the upper arm (Fig. [7] bottom right). 
Such structures can thus be used to rewire the skyrmion 
current into specific parts of the circuit. 


VI. SUMMARY 

It is shown that an electric-field provides an addi¬ 
tional degree of freedom for effective control of current- 
driven skyrmion motion. By using a symmetry-based 



















phenomenology we find the general form of spin-orbit- 
induced electric-field coupling for thin films with broken 
inversion symmetry about the film plane. In general, 
this coupling can be classified as dissipative, i.e. not 
derivable from a free energy, and conservative, i.e. terms 
which can be understood as electric-field- induced mod¬ 
ification of the free energy. The conservative terms can 
further be identified as electric-field dependent perpen¬ 
dicular anisotropy (when electric field does not couple to 
gradients of magnetization) and electric-field dependent 
exchange and Dzyaloshinskii-Moriya interaction (when 
electric-field couples to gradient in magnetization). Mo¬ 
tivated by applications and available experimental data, 
we specifically studied effects of the electric-field depen¬ 
dent perpendicular ansisotropy in thin films, which also 
harbor a nanosized skyrmion state stabilized by dipo¬ 
lar interaction. We find that there exist two regimes 
for controlling the static and hence the dynamical prop¬ 
erties of skyrmions: (a) below a critical value, such 
that the corresponding perpendicular ansitoropy is still 
strong enough to overcome the out-of-plane dipolar in¬ 
teraction, skyrmions are a stable configuration showing a 
decrease in diameter and domain wall width surrounding 
the skyrmion with increasing anisotropy, and (b) above 
the critical value, the skyrmion itself becomes unstable 
due to formation of multi-domain/ in-plane magnetized 
state. The decrease in diameter and domain wall width 
with increasing anisotropy in regime (a) is, in turn, driven 
by an attempt to reduce the wall energy. This change 
of static properties modulates both the speed and the 
angle of skyrmion motion driven by a current induced 
spin-orbit force. This modulation can be understood 
within Thiele’s analytical model of rigid skyrmion mo¬ 
tion, where, the static properties enter the velocity via 
D s and A w dependence of spin-orbit and damping forces 


acting on the skyrmion. On the other hand, correspond¬ 
ing to regime (b), desired unstable regions can be defined 
for skyrmions using gate voltages, which are observed 
in micromagnetic simulations to guide skyrmions along 
the boundary between stable and unstable regions.This 
regime of skyrmion guiding is established by appearance 
of additional forces due to variation of skyrmion’s poten¬ 
tial energy, caused by the dipolar interactions, as the 
skyrmion approaches the unstable region. Using this 
ability to electrically design potential energy landscapes, 
we demonstrate transistor-like and multiplexer-like spin- 
tronic devices controlling skyrmion current by defining 
electric-field gates. The proposed electric-field control is 
thus expected to extend the reach of skyrmionics into 
logic and computing devices. 
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